Method for compensating a magnetic locator, locator and computer program

ABSTRACT

The invention relates to a method for compensating a magnetic locator in the presence of a magnetic-field-disturbing material, comprising: an emitter ( 10 ) comprising at least one coil emitting an emission magnetic field; a receiver ( 20 ) comprising at least one receiving coil and a device providing a plurality of measurements Ip i  of a receiving magnetic field induced by the emission field in each receiving coil; and a processing unit ( 25 ) comprising a field model allowing the calculation of a position (P) and/or an orientation (Q) of the receiver by means of calculation of a prediction H i  of the measurements according to a criterion (C) calculated according to an error E i  which is itself calculated in relation to the measurements Ip i . The invention is characterised in that the error E i  is calculated by successive iterations from initial values prescribed by the prediction H i  as being the difference between the measurements Ip i  and a disturbed model Hp i , according to the equation E i =Ip i −Hp i , the disturbed model Hp i  satisfying Hp i =H i +P i  (α i =−arctan(βω i ), (I) the parameter β being identical for all of the measurements Ip i , the calculation being carried out in such a way as to minimise the criterion C.

The invention relates to a method for compensating a magnetic locator in the presence of at least one disruptive material.

A magnetic locator generally comprises a transmitting device, the so-called transmitter, having one or more transmitter coils rigidly connected to one another, and a receiving device, the so-called receiver, having one or more receiver coils rigidly connected to one another. The joint analysis of the magnetic fields transmitted by the transmitter coils and of the magnetic fields measured by the receiver coils makes it possible to determine by a processing unit the position and/or orientation of the receiving device with respect to the transmitting device. This processing unit comprises a field model making it possible to compute the position and/or orientation of the receiver by computing a prediction of the measurements as a function of a criteria of minimization of an error computed with respect to the measurements.

The problem of the invention is a material can exist that disrupts the magnetic field between the transmitter and the receiver. In particular, certain electrically conductive materials can give rise to eddy currents when these materials are placed in a magnetic field. The eddy currents that then circulate in this conductive material in turn generate a disruptive magnetic field.

Thus, the presence of a disruptive material causes errors in the computation of the prediction by the processing unit.

The invention aims to solve this problem by proposing a method and a device for compensating a magnetic locator in the presence of at least one disruptive material, which make it possible to compute the correct position and/or orientation of the receiver corresponding as much as possible to its actual position and/or to its actual orientation, despite the disruptive magnetic field induced by the material.

For this purpose, a first subject of the invention is a method for compensating a magnetic locator in the presence of at least one magnetic-field-disrupting material, said magnetic locator comprising:

-   -   at least one transmitter comprising at least one transmitter         coil and at least one generator of at least one transmission         signal, connected to the at least one transmitter coil, so that         the at least one transmitter coil transmits at least one         transmitting magnetic field at at least one determined frequency         ω_(i) in response to the transmission signal that is sent to it         by the generator,     -   at least one receiver comprising at least one receiver coil and         a measuring device, which is connected to the at least one         receiver coil and which supplies at least one measurement Ip_(i)         of a receiving magnetic field respectively induced by the         transmitting magnetic field in each receiver coil, in such a way         as to supply several measurements Ip_(i) for i ranging from 1 to         N,     -   a processing unit comprising a field model making it possible to         compute a position and/or an orientation of the receiver by         computing a prediction H_(i) of the measurements as a function         of a criterion computed as a function of an error E_(i), itself         computed with respect to the measurements Ip_(i),

characterized in that

the error E_(i) is computed by successive iterations from initial prescribed values of the prediction H_(i) as being the difference between the measurements Ip_(i) and a disrupted model Hp_(i) of the measurements, according to the equation E _(i) =Ip _(i) −Hp _(i),

the disrupted model Hp_(i) of the measurements verifying the following equations

Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢ e j ⁢ α i ⁢ α i = - arctan ⁡ ( β · ω i ) ,

where P_(i) is the disruption made to the measurements Ip_(i) by the magnetic-field-disrupting material, ρ_(i) is the intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, β being a parameter of the magnetic-field-disrupting material,

the parameter β being identical for all the measurements Ip_(i),

the computation being carried out in such a way as to minimize the criterion.

According to an embodiment of the invention, the field model used to compute a position and/or an orientation of the receiver by computing the prediction H_(i) of the measurements as a function of the criterion uses a Levenberg-Marquardt minimization algorithm.

A number of first embodiments are described below. For each embodiment, the following steps described below are carried out in each iteration.

According to a first embodiment of the invention, at each iteration,

-   -   the prediction H_(i) is initialized at the initial prescribed         values,     -   then are computed Δ_(i), the parameter β so as to minimize the         vector A·β−B over its coordinates and the disrupted model Hp_(i)         as a function of H_(i) and Ip_(i) according to the equations         Δ_(i)=[         (H _(i))·         (Ip _(i))−         (H _(i))·         (Ip _(i))]/[         (H _(i))²+         (H _(i))²]         Hp _(i) =H _(i)·(1+Δ_(i)·(j+βω _(i))),

where A is a vector, the first and second coordinates of which are respectively formed by: A_(2i)=Δ_(i)·

(H_(i))·ω_(i) et A_(2i+1)=Δ_(i)·

(H_(i))·ω_(i),

B is a vector, the first and second coordinates of which are respectively formed by: B _(2i)=

(Ip _(i))−

(H _(i))+Δ_(i)·

(H _(i)) and B _(2i+1)=

(Ip _(i))−

(H _(i))−Δ_(i)·

(H _(i)),

for i ranging from 1 to N, where

(H_(i)) is the real part of the prediction H_(i),

(H_(i)) is the imaginary part of the prediction H_(i),

(Ip_(i)) is the real part of the measurement Ip_(i),

(Ip_(i)) is the imaginary part of the measurement Ip_(i),

-   -   then the error E_(i) is computed,     -   then the criterion C is computed according to the invention

${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$

-   -   then the prediction H_(i) corresponding to the criterion C is         computed by the field model,

this computed prediction H_(i) being used for the following iteration until the criterion C becomes less than a prescribed non-zero positive bound η.

Following another first embodiment of the invention, at each iteration

-   -   the prediction H_(i) is initialized at the prescribed initial         values,     -   then are computed Δ_(i), the parameter β and the disrupted model         Hp_(i) as a function of H_(i) and Ip_(i) according to the         equations

${{\Delta_{i} = {\left\lbrack {{{\left( H_{i} \right) \cdot}\left( {Ip}_{i} \right)} - {{\left( H_{i} \right) \cdot}\left( {Ip}_{i} \right)}} \right\rbrack/\left\lbrack \left( {\left( H_{i} \right)^{2} + \left( H_{i} \right)^{2}} \right. \right\rbrack}},{\beta = \frac{A^{T} \cdot B}{{❘A❘}^{2}}}}{{{Hp}_{i} = {H_{i} \cdot \left( {1 + {\Delta_{i} \cdot \left( {j + {\beta \cdot \omega_{i}}} \right)}} \right)}},}$

where A is a vector, the first and second coordinates of which are respectively formed by: A_(2i)=Δ_(i)·

(H_(i))·ω_(i) et A_(2i+1)=Δ_(i)·

(H_(i))·ω_(i),

B is a vector, the first and second coordinates of which are respectively formed by:

B_(2i)=

(Ip_(i))−

(H_(i))+Δ_(i)·

(H_(i)) and B_(2i+1)=

(Ip_(i))−

(H_(i))−Δ_(i)·

(H_(i)), for i ranging from 1 to N, where

(H_(i)) is the real part of the prediction H_(i),

(H_(i)) is the imaginary part of the prediction H_(i),

(Ip_(i)) is the real part of the measurement Ip_(i),

(Ip_(i)) is the imaginary part of the measurement Ip_(i).

-   -   then the error E_(i) is computed,     -   then the criterion C is computed according to the equation

${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$

-   -   then the field model is used to compute the prediction H_(i)         corresponding to the criterion C,

this computed prediction H_(i) being used for the following iteration until the computed criterion C becomes less than a prescribed non-zero positive bound η.

A number of second embodiments are described below.

According to a second embodiment of the invention, at each iteration

-   -   the prediction H_(i) is initialized at the initial prescribed         values,     -   then the parameter β is computed so as to minimize the vector         A′·β+B′ on its coordinates, then the Δ_(i) and the disrupted         model Hp_(i) as a function of H_(i) and Ip_(i) according to the         equations:

${\Delta_{i} = \frac{\left( {Ip}_{i} \right) - \left( H_{i} \right)}{\left( H_{i} \right) + {\left( H_{i} \right)\beta\omega_{i}}}}{{Hp}_{i} = {H_{i} \cdot \left( {1 + {\Delta_{i} \cdot \left( {j + {\beta \cdot \omega_{i}}} \right)}} \right)}}$

where A′ is a vector, the coordinates of which are respectively formed by: ω_(i)·[

(H_(i))·

(Ip_(i))−

(Ip_(i))·

(H_(i))],

B′ is a vector, the coordinates of which are respectively formed by:

(H_(i))²+

(H_(i))²−(

(H_(i))·

(Ip_(i))+

(H_(i))·(Ip_(i))),

for i ranging from 1 to N, where

(H_(i)) is the real part of the prediction H_(i),

(H_(i)) is the imaginary part of the prediction H_(i),

(Ip_(i)) is the real part of the measurement Ip_(i),

(Ip_(i)) is the imaginary part of the measurement Ip_(i),

-   -   then the error E_(i) is computed,     -   then the criterion C is computed according to the equation

${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$

-   -   then the field model is used to compute the prediction H_(i)         corresponding to the criterion C,

this computed prediction H_(i) being used for the following iteration until the computed criterion C becomes less than a prescribed non-zero positive bound η.

According to another second embodiment of the invention, at each iteration

-   -   the prediction H_(i) is initialized at the initial prescribed         values,     -   then the parameter β is computed, then the Δ_(i) and the         disrupted model Hp_(i) as a function of H_(i) and Ip_(i)         according to the equations:

${\beta = \frac{A^{\prime T} \cdot B^{\prime}}{{❘A^{\prime}❘}^{2}}}{\Delta_{i} = \frac{\left( {Ip}_{i} \right) - \left( H_{i} \right)}{\left( H_{i} \right) + {\left( H_{i} \right)\beta\omega_{i}}}}{{Hp}_{i} = {H_{i} \cdot \left( {1 + {\Delta_{i} \cdot \left( {j + {\beta \cdot \omega_{i}}} \right)}} \right)}}$

where A′ is a vector, the coordinates of which are respectively formed by: ω_(i)·[

(H_(i))·

(Ip_(i))−

(Ip_(i))·

(H_(i))],

B′ is a vector, the coordinates of which are respectively formed by:

(H_(i))²+

(H_(i))²−(

(H_(i))·

(Ip_(i))+

(H_(i))·(Ip_(i))),

for i ranging from 1 to N, where

(H_(i)) is the real part of the prediction H_(i),

(H_(i)) is the imaginary part of the prediction H_(i),

(Ip_(i)) is the real part of the measurement Ip_(i),

(Ip_(i)) is the imaginary part of the measurement Ip_(i),

-   -   then the error E_(i) is computed,     -   then the criterion C is computed according to the equation

${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$

-   -   then the field model is used to compute the prediction H_(i)         corresponding to the criterion C,

this computed prediction H_(i) being used for the following iteration until the computed criterion C becomes less than a prescribed non-zero positive bound η.

A number of third embodiments are described below.

According to a third embodiment of the invention, at each iteration

-   -   the prediction H_(i) is initialized at the prescribed initial         values,     -   the parameter β is then added as a state variable of the         Levenberg-Marquardt minimization algorithm,     -   then the Δ_(i) are computed as a function of H_(i) and Ip_(i)         and β according to the equations:         Hp _(i) =H _(i)·(1+Δ_(i)·(j+βω _(i))),         the Δ_(i) being a solution to the following system of equations:         (H _(i))−         (Ip _(i))+Δ_(i)·(         (H _(i))·β·ω_(i)−         (H _(i)))=0         and         (H _(i))−         (Ip _(i))+Δ_(i)·(         (H _(i))+         (H _(i))·β·ω_(i))=0         for i ranging from 1 to N, where

(H_(i)) is the real part of the prediction H_(i),

(H_(i)) is the imaginary part of the prediction H_(i),

(Ip_(i)) is the real part of the measurement Ip_(i),

(Ip_(i)) is the imaginary part of the measurement Ip_(i),

-   -   then the error E_(i) is computed,     -   then the criterion C is computed according to the equation

${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$

-   -   then the field model is used to compute the prediction H_(i)         corresponding to the criterion C,

this computed prediction H_(i) being used for the following iteration until the computed criterion C becomes less than a prescribed non-zero positive bound η.

According to a third embodiment of the invention, at each iteration

-   -   the prediction H_(i) is initialized at the prescribed initial         values,     -   the parameter β is then added as a state variable of the         Levenberg-Marquardt minimization algorithm,     -   then the Δ_(i) are computed as minimizing a vector C·Δ-D over         its coordinates as a function of H_(i) and Ip_(i) and β         according to the equations:         Hp _(i) =H _(i)·(1+Δ_(i)·(j+β·ω _(i))),

where Δ is a vector Δ having Δ_(i) as coordinates,

D is a vector, the first and second coordinates of which are respectively formed by: −(

(H_(i))−

(Ip_(i))) and −(

(H_(i))−

(Ip_(i))),

C is a matrix, having as coefficients corresponding to Δ_(i) respectively

(H_(i))·β·ω_(i)−

(H_(i)) and

(H_(i))+

(H_(i))·β·ω_(i) and the coefficients 0 elsewhere, for i ranging from 1 to N, where

(H_(i)) is the real part of the prediction H_(i),

(H_(i)) is the imaginary part of the prediction H_(i),

(Ip_(i)) is the real part of the measurement Ip_(i),

(Ip_(i)) is the imaginary part of the measurement Ip_(i),

-   -   then the error E_(i) is computed,     -   then the criterion C is computed according to the equation

${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$

-   -   then the field model is used to compute the prediction H_(i)         corresponding to the criterion C,         this computed prediction H_(i) being used for the following         iteration until the computed criterion C becomes less than a         prescribed non-zero positive bound η.

According to an embodiment of the invention, the determined frequencies ω_(i) are separate from one another.

According to an embodiment of the invention, the transmitter comprises a plurality of transmitter coils respectively transmitting a plurality of transmitting magnetic fields,

the measuring device measuring for each receiver coil a plurality of receiving magnetic fields respectively induced by the plurality of transmitting magnetic fields in the receiver coil and forming the measurements Ip_(i) for i ranging from 1 to N.

According to another embodiment of the invention, the locator comprises a transmitter coil transmitting at K different frequencies ω_(jk) and J receiver coils, the measuring device supplying the measurements Ip_(jk) for the index j ranging from 1 to J and the index k ranging from 1 to K,

the processing unit computes

ρ j = ∑ k = 1 K u jk ⁣ · ( ❘ "\[LeftBracketingBar]" H jk ❘ "\[RightBracketingBar]" 2 + ❘ "\[LeftBracketingBar]" Ip jk ❘ "\[RightBracketingBar]" 2 - 2 · v jk ) ∑ k = 1 K u jk 2 ⁢ and ⁢ β = ∑ j = 1 J ∑ k = 1 K u jk · ( v jk - ❘ "\[LeftBracketingBar]" H jk ❘ "\[RightBracketingBar]" 2 ) ∑ j = 1 J ∑ k = 1 K u jk 2 ⁢ with ⁢ v jk = ( H jk ) · ⁢ ( Ip jk ) + ( H jk ) · ⁢ ( Ip jk ) ⁢ u jk = ( ( H jk ) · ⁢ ( Ip jk ) - ( H jk ) · ⁢ ( Ip jk ) ) · ω jk ⁢ Ip jk = H jk · ( 1 + Δ jk · ( j + β · ω jk ) ) ⁢ Δ jk = ρ j ⁢ ω j ⁢ k 1 + ( β · ω jk ) 2

where H_(jk) is the prediction, ρ_(j) is the intensity of the disruption, Hp_(jk) is the disrupted model and Hp_(jk)−Ip_(jk)=0.

According to another embodiment of the invention, the locator comprises L transmitter coils or L transmitters, which transmit at K different frequencies ω_(jkl), and J receiver coils, the measuring device supplying the measurements Ip_(jkl) for the index j ranging from 1 to J, the index k ranging from 1 to K and the index l ranging from 1 to L,

the processing unit computes

${\rho_{jl} = \frac{\sum\limits_{k = 1}^{K}{u_{jkl} \cdot \left( {{❘H_{jkl}❘}^{2} + {❘{Ip}_{jkl}❘}^{2} - {2 \cdot v_{jkl}}} \right)}}{\sum\limits_{k = 1}^{K}u_{jkl}^{2}}}{and}{\beta_{l} = \frac{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}{u_{jkl} \cdot \left( {v_{jkl} - {❘H_{jkl}❘}^{2}} \right)}}}{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}u_{jkl}^{2}}}}{with}{v_{jkl} = {{{\left( H_{jkl} \right) \cdot}\left( {Ip}_{jkl} \right)} + {{\mathcal{J}\left( H_{jkl} \right)} \cdot {\mathcal{J}\left( {Ip}_{jkl} \right)}}}}{u_{jkl} = {\left( {{\left( H_{jkl} \right) \cdot {\mathcal{J}\left( {Ip}_{jkl} \right)}} - {{{\mathcal{J}\left( H_{jkl} \right)} \cdot}\left( {Ip}_{jkl} \right)}} \right) \cdot \omega_{jkl}}}{{Ip}_{jkl} = {H_{jkl} \cdot \left( {1 + {\Delta_{jkl} \cdot \left( {j + {\beta \cdot \omega_{jkl}}} \right)}} \right)}}{\Delta_{jkl} = {\rho_{j}\frac{\omega_{jkl}}{1 + \left( {\beta \cdot \omega_{jkl}} \right)^{2}}}}$

where H_(jkl) is the prediction, ρ_(jl) is the intensity of the disruption, Hp_(jkl) is the disrupted model and Hp_(jkl)−Ip_(jkl)=0.

A second subject of the invention is a magnetic locator comprising:

-   -   at least one transmitter comprising at least one transmitter         coil and at least one generator able to generate at least one         transmission signal, connected to the at least one transmitter         coil, so that the at least one transmitter coil is able to         transmit at least one transmitting magnetic field at at least         one determined frequency ω_(i) in response to the transmission         signal that is sent to it by the generator,     -   at least one receiver comprising at least one receiver coil and         a measuring device, which is able to be connected to the at         least one receiver coil and which is able to supply at least one         measurement Ip_(i) of a receiving magnetic field respectively         induced by the transmitting magnetic field in each receiver         coil, to supply several measurements Ip_(i) for i ranging from 1         to N,     -   a processing unit comprising a field model making it possible to         compute a position and/or an orientation of the receiver by         computing a prediction H_(i) of the measurements as a function         of a criterion computed as a function of an error E_(i), itself         computed with respect to the measurements Ip_(i),

characterized in that

the processing unit is configured so that the error E_(i) is computed by successive iterations from prescribed initial values of the prediction H_(i) as being the difference between the measurements Ip_(i) and a disrupted model Hp_(i) of the measurements, according to the equation E _(i) =Ip _(i) −Hp _(i),

the disrupted model Hp_(i) of the measurements verifying the following equations

Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢ e jα i , α i = - arctan ⁡ ( β · ω i ) ,

where P_(i) is a disruption made to the measurements Ip_(i) by a magnetic-field-disrupting material, ρ_(i) is the intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, β being a parameter of the magnetic-field-disrupting material,

the parameter β being identical for all the measurements Ip_(i),

the processing unit being configured to carry out the computation in such a way as to minimize the criterion.

The locator may be produced according to several embodiments, where the processing unit is configured to implement the steps of respectively each embodiment described above of the method at each iteration.

A third subject of the invention is a computer program, comprising code instructions for implementing the method for compensating a magnetic locator in the presence of at least one magnetic-field-disrupting material, as described above, when the computer program is executed on a calculator.

The invention will be better understood on reading the following description, given solely by way of non-limiting example with reference to the appended drawings, wherein:

FIG. 1 schematically represents a modular block diagram of a magnetic locator according to an embodiment of the invention,

FIG. 2 schematically represents a magnetic-field-disrupting material,

FIGS. 3, 4 and 5 schematically represent three magnetic phase diagrams for three different materials,

FIGS. 6 and 7 are Fresnel diagrams representing the disrupted measurement and a disrupted model Hp_(i), used according to an embodiment of the invention,

FIG. 8 is a modular block diagram of a processing unit of the locator according to an embodiment of the invention,

FIG. 9 is a modular block diagram of a processing unit of the locator according to FIG. 8, using a first algorithm,

FIG. 10 is a modular block diagram of a processing unit of the locator according to FIG. 8, using a second algorithm,

FIG. 11 is a modular block diagram of a processing unit of the locator according to FIG. 8, using a third algorithm.

Below is a description, with reference to the figures, of a method for compensating a magnetic locator and a device for compensating a magnetic locator.

In FIG. 1, a magnetic locator 1 or device 1 for magnetic location comprises one or more transmitters 10 and one or more receivers 20.

The transmitter 10 comprises a transmitter coil or a plurality (number Ne) of transmitter coils, for example E1, E2, E3, and a generator or a plurality of generators, for example 101, 102, 103, of at least one transmission signal, connected to the transmitter coils E1, E2, E3, to control the transmitting coils with their respective transmission signal. The transmitter coils are rigidly attached to one and the same mechanical support of the transmitter 10, in order to occupy a known position and have a known orientation and can for example be oriented along separate axes. The transmitter 10 has for example a known position and a known orientation. Here, by way of non-limiting example, the case is considered of a system in which the transmitter 10 comprises 3 transmitter coils E1, E2, E3 oriented along three separate axes, for example along the three axes of an orthonormal reference frame. These three axes of the respective magnetic moments of the coils E1, E2, E3 can be colinear with the axis of the winding of the coils in the case of a coil orthogonal to the axis of winding. They can however be different therefrom if the coil is not at 90°, for example at 45°. When a voltage is imposed across the terminals of a transmitter coil E1, E2, E3, a current circulates in the transmitter coil E1, E2, E3 which then generates a magnetic field Be1, Be2, Be3 proportional to the current through it and the shape of which depends on the features of the coil (orientation, magnetic moment, shape etc.). A voltage source may, for example, be used as a generator 101, 102, 103 to impose a voltage across the terminals of the transmitter coil E1, E2, E3 driving the creation of a current.

When each generator 101, 102, 103 is started and connected to the transmitter coil E1, E2, E3 respectively, each transmitter coil E1, E2, E3 transmits a transmitting magnetic field, Be1, Be2, Be3 respectively, in response to the transmission signal that is sent to it by the generator 101, 102, 103. The transmission signal and the transmitting magnetic field, Be1, Be2, Be3 respectively, are for example sinusoidal and have for example a determined frequency, f1, f2, f3 respectively. The determined frequencies may be identical to one another or different from one another. The determined frequencies may for example be of 100 Hz to 50 kHz.

The receiver 20 comprises a receiver coil or several (number Nr) receiver coils, for example R1, R2, R3, and a measuring device 21, which is connected to the receiver coils R1, R2, R3. The number Nr of receiver coils may be different from or equal to the number Ne of transmitter coils. The receiver coils are rigidly attached to one and the same mechanical support of the receiver 120 and may for example be oriented along separate axes. The receiver 20 has for example a known position and a known orientation. Here, by way of example, the case is considered of a system in which the receiver 20 comprises 3 receiver coils R1, R2, R3 oriented along three separate axes, for example along the three axes of an orthonormal frame of reference. These three axes of the respective magnetic moments of the coils R1, R2, R3 may be colinear with the axis of the winding of the coils in the case of a coil orthogonal to the axis of winding. They may however be different therefrom if the coil is not at 90°, for example at 45°. In the presence of a variable magnetic field, a voltage (measurement) proportional to the variation of the flux of the magnetic field appears in the receiver coil. By measuring the voltage across the terminals of the receiver coil, for example using a voltmeter or other similar means, or by measuring the current (measurement) flowing through the receiver coil, for example by an ampermeter or other similar means, it is possible to determine the magnetic field to which the receiver coil is subjected, on the condition that the features of the receiver coil, which particularly comprise the magnetic moment of the receiver coil, are known.

The measuring device 24 measures for each receiver coil R1, R2, R3 a plurality of receiving magnetic fields respectively induced by the plurality of transmitting magnetic fields Be1, Be2, Be3 and called measurements Ip_(i) for i ranging from 1 to N and comprises for example for this purpose measuring modules 21, 22, 23 respectively connected to the receiver coils R1, R2, R3. We have N=Nr·Ne possible pairs of transmitter coils and receiver coils, which makes it possible to constitute N measurements Ip_(i) by the measuring device 24. According to an embodiment, N is greater than 1. The measurements Ip_(i) are for example complex, each having a gain and a phase. The measurements Ip_(i) are for example supplied on a demodulation output of the measuring device 24. In the case where the transmission signal is sinusoidal with a determined frequency, the corresponding measurement Ip_(i) is also sinusoidal with the same determined frequency. Thus, in the example of FIG. 1, the measuring device 21 measures via the measuring module 21 connected to the receiver coil R1:

-   -   the receiving magnetic field Br11 (written E1→R1), which is         induced in the receiver coil R1 by the transmitting magnetic         field Be1 of the transmitter coil E1,     -   the receiving magnetic field Br21 (written E2→R1), which is         induced in the receiver coil R1 by the transmitting magnetic         field Be2 of the transmitter coil E2,     -   the receiving magnetic field Br31 (written E3→R1), which is         induced in the receiver coil R1 by the transmitting magnetic         field Be3 of the transmitting coil E1.

The measuring device 21 measures via the measuring module 23 connected to the receiver coil R2:

-   -   the receiving magnetic field Br12 (written E1→R2), which is         induced in the receiver coil R2 by the transmitting magnetic         field Be1 of the transmitter coil E1,     -   the receiving magnetic field Br22 (written E2→R2), which is         induced in the receiver coil R2 by the transmitting magnetic         field Be2 of the transmitter coil E2,     -   the receiving magnetic field Br32 (written E3→R2), which is         induced in the receiver coil R2 by the transmitting magnetic         field Be3 of the transmitter coil E1.

The measuring device 21 measures for the measuring module 22 connected to the receiver coil R3:

-   -   the receiving magnetic field Br13 (written E1→R3), which is         induced in the receiver coil R3 by the transmitting magnetic         field Be1 of the transmitter coil E1,     -   the receiving magnetic field Br23 (written E2→R3), which is         induced in the receiver coil R3 by the transmitting magnetic         field Be2 of the transmitter coil E2,     -   the receiving magnetic field Br33 (denoted E3→R3), which is         induced in the receiver coil R3 by the transmitting magnetic         field Be3 of the transmitter coil E1.

The transmitter 10 and the receiver 20 are connected to a processing unit 25, (which can for example comprise a microprocessor) connected on the one hand to the generator 101, 102, 103 of the transmitting device 10 and on the other hand to the measuring device 24 of the receiving device 20. The processing unit 25 is configured to process the transmission signals transmitted and the measurements Ip_(i) received, and thus makes it possible to retrieve from these transmitted transmission signals and from these measurements Ip_(i) the position P and the orientation Q of the receiver 20 relative to the transmitter 10. The processing unit 25 receives information concerning the features of the transmission signals applied to the transmitter coils as well as information concerning the features of the signals flowing through the receiver coils (representative of the receiving magnetic field by the receiver coil and therefore the measurements Ip_(i)). Based on the intensity and the phase of the receiving magnetic fields sensed by the receiver coils, the processing unit 25 determines the spatial coordinates of position P and orientation Q of all the receiver coils with respect to all the transmitter coils, for example using a minimization algorithm, for example of Levenberg-Marquardt type, making it possible to minimize the error between the measured fields (receiving magnetic fields or measurements Ip_(i)) and the theoretical magnetic fields modelled on the basis of a priori knowledge about the transmitter and the receiver.

A magnetic-field-disrupting material 3 can be present near the locator 1, for example between the transmitter 10 and the receiver 20 or near them, as represented in FIG. 2, where only a transmitter coil of the transmitter 10 and a receiver coil of the receiver 20 have been represented by way of example. The disruptive material 3 may be for example an electrically conductive material, for example metal. Any electrically conductive material 3 placed in a variable magnetic field gives rise to eddy currents in the conductive material 3, symbolized by the current loops CF in FIG. 2. The eddy currents CF that then circulate in this conductor 3 in turn give rise to a disruptive magnetic field, the phase of which depends on the impedance of the conductor 3. But this phase is constant for a given frequency. The disruptive material 3 can for example be a metal prosthesis, such as for example a femoral head of a hip prosthesis, a femoral stem of a hip prosthesis, a metal plate, for example made of copper, another metal coil (other than the transmitter coils and the receiver coils) or other.

If the magnetic field induced by the disruptor is not in phase with the signal of the transmitter then it is possible to detect the disruption and to partly compensate for it. The outlined method makes it possible to compensate for a single disruptor, it does not compensate for disruptions in phase with the signal.

The situations can be summarized by the three phase diagrams of FIGS. 3, 4 and 5 (real part Re of the magnetic fields on the abscissae, imaginary part Im of the magnetic fields on the ordinates), with Bp the initial magnetic field without disruption (for example the transmitting magnetic field Be1 or Be2 or Be3), Bs the magnetic field transmitted by the eddy currents CF of the magnetic-field-disrupting material 3 and Bc the resulting magnetic field, equal to the sum of Bp and Bs.

In FIG. 3, in the case of a ferromagnetic material (such as for example ferrites), the material generates a ferromagnetic disruption, where there are no eddy currents and there is no phase shift between the initial field Bp and the resulting field Bc.

In FIG. 4, in the case of a magnetic-field-disrupting material 3, which is electrically conductive (such as for example a copper and/or silver and/or aluminum plate), the eddy currents CF induced by the material 3 create a secondary field Bs which has the effect of reducing the initial field Bp and introducing a phase shift. The resulting field Bc is phase-shifted and lower in amplitude with respect to the initial field Bp.

In FIG. 5, in the case of a magnetic-field-disrupting material 3, which is mixed (both ferromagnetic and electrically conductive), for relatively low frequencies, the ferromagnetic nature means that the amplitude of the resulting field Bc is higher than that of the initial field Bp. The eddy currents CF induced by the material 3 at the same time introduce a phase shift. As the frequency increases, the effect of the ferromagnetism decreases, that of the eddy currents CF induced by the material 3 CF increases, the amplitude of the resulting field Bc decreases more and more, the phase shift increases then decreases according to a law specific to the material 3.

The invention makes it possible to compensate for the disruption Bs due to the eddy currents CF induced by the material 3. The ferromagnetic disruptions are assumed to be zero.

The magnetic-field-disrupting material 3 can be a conductive metal object at a certain position in space in relation to a transmitter coil. This implies that if the metal object moves, it constitutes a new disruptor, the features of which are different.

The magnetic locator 1 implements the method for compensating for the presence of the magnetic-field-disrupting material 3, described below.

In FIG. 8, the locator 1 comprises the processing unit 25 comprising a module 26 for modelling the field making it possible to compute a position P and/or an orientation Q of the receiver 20 by computing a prediction H_(i) of the measurements as a function of a criterion C of minimization of an error E_(i) computed with respect to the measurements Ip_(i). This module 26 may for example comprise a Kalman filter for computing the prediction H_(i). The module 26 disposes of a physical model of the locator 1, which predicts in the absence of disruption and as a function of the position P and/or the orientation Q of the receiver 20 the N measurements, to supply the predictions H_(i). The module 26 searches for the position P and/or the orientation Q (or state of the sensor P, Q) and computes the prediction H_(i) verifying the criterion C. The criterion C is considered by the module 26 as representing the deviation (which can for example be a norm of the error taken in the meaning of L2, which can be the sum of the squares of the errors E_(i) over i ranging from 1 to N) between the prediction H_(i) and the measurements Ip_(i). The criterion C can be computed as being for example a norm of the error E_(i) taken in the meaning of L². According to an embodiment, the criterion C can be computed as being the sum of the squares of the errors E_(i) over i ranging from 1 to N. According to an embodiment, the criterion C can be computed according to the following equation

$C = {\sum\limits_{i = 0}^{N}{❘{{H_{i}\left( {P,Q} \right)} - {IP}_{i}}❘}^{2}}$

According to an embodiment, the error E_(i) is computed by an error computing module 27 as being the difference between the measurements Ip_(i) and a disrupted model Hp_(i) of the measurements, according to the equation E _(i) =Ip _(i) −Hp _(i).

At a given instant, there are N complex measurements Ip_(i) and N complex predictions H_(i) resulting from the model 26, for i ranging from 1 to N. Each of these measurements Ip_(i) is made at a certain frequency or angular frequency (equal to the frequency multiplied by 2π), which is written below ω_(i).

The measurements Ip_(i) are for example stored in a vector Ip having as coordinates these measurements Ip_(i), i.e.

${Ip} = {\begin{pmatrix} {{Br}11} \\ {{Br}21} \\ \ldots \\ \frac{{BrN}_{r}1}{{Br}21} \\ \ldots \\ \ldots \\ {{BrN}_{e}N_{r}} \end{pmatrix} = \begin{pmatrix} \left. {E1}\rightarrow{R1} \right. \\ \left. {E1}\rightarrow{R2} \right. \\ \ldots \\ \frac{\left. {E1}\rightarrow R_{N_{r}} \right.}{\left. {E2}\rightarrow{R1} \right.} \\ \ldots \\ \ldots \\ \left. E_{N_{e}}\rightarrow R_{N_{r}} \right. \end{pmatrix}}$

The predictions H_(i) are for example stored in a vector H(P,Q) having as coordinates these predictions H_(i), i.e.

${H\left( {P,Q} \right)} = {\begin{pmatrix} {H\left( {{Br}11} \right)} \\ {H\left( {{Br}21} \right)} \\ \ldots \\ \frac{H\left( {{BrN}_{r}1} \right)}{H\left( {{Br}21} \right)} \\ \ldots \\ \ldots \\ {H\left( {{BrN}_{e}N_{r}} \right)} \end{pmatrix} = \begin{pmatrix} {H\left( {E1}\rightarrow{R1} \right)} \\ {H\left( {E1}\rightarrow{R2} \right)} \\ \ldots \\ \frac{H\left( {E1}\rightarrow R_{N_{r}} \right)}{H\left( {E2}\rightarrow{R1} \right)} \\ \ldots \\ \ldots \\ {H\left( E_{N_{e}}\rightarrow R_{N_{r}} \right)} \end{pmatrix}}$

In the absence of any magnetic-field-disrupting material 3 by eddy current P_(i)=0, it is supposed that the model 26 is appropriately calibrated and that it is correctly following the measurement (then called I_(i) in the absence of any magnetic-field-disrupting material 3) as a function of the time t, according to the equation H_(i)(t)=I_(i)(t)+ε_(i) with ε_(i) negligible (ε_(i)=0 in the remainder of the text).

When a magnetic-field-disrupting material 3 is present, a disruption P_(i), non-negligible, is added to the measurement I_(i) and this then gives Ip_(i)=I_(i)+P_(i) and the disrupted model Hp_(i) is equal to Hp_(i)=H_(i)+P_(i). FIG. 6 is a Fresnel diagram representing the disrupted measurement Ip_(i) with respect to a non-disrupted measurement I_(i). FIG. 7 is a Fresnel diagram representing the non-disrupted prediction H_(i) with respect to the disrupted model Hp_(i).

No correction is made to the measurement Ip_(i) which is itself accurate, but the disruption P_(i) is modeled by computing a disrupted model Hp_(i) such that the criterion C is minimal. According to an embodiment, the processing unit 25 comprises a module 29 for computing the criterion C of minimization of an error E_(i) computed with respect to the measurements Ip_(i). According to an embodiment,

$C = {\sum\limits_{i = 0}^{N}{❘{{HP}_{i} - {IP}_{i}}❘}^{2}}$

According to an embodiment, one computes as disrupted model Hp_(i) a projection of the prediction H_(i) onto the measurements Ip_(i) as a function of a parameter β characteristic of the magnetic-field-disrupting material 3, and for example only as a function of the parameter β. This projection makes it possible to reduce the effect of the magnetic-field-disrupting material 3 by eddy currents.

According to an embodiment, the criterion C=E^(T)·E is minimized with E=Ip−Hp(H,Ip,β).

The vector Hp(H,Ip,β) is a projection function of H onto Ip and is computed as a function of the vector H(P, Q) and of the vector Ip.

The vector Hp(H,Ip,β) has as coordinates this disrupted model Hp_(i), for i ranging from 1 to N, i.e.

${{Hp}\left( {H,{Ip},\beta} \right)} = {\begin{pmatrix} {{Hp}\left( {{Br}11} \right)} \\ {{Hp}\left( {{Br}21} \right)} \\ \ldots \\ \frac{{Hp}\left( {{BrN}_{r}1} \right)}{{Hp}\left( {{Br}21} \right)} \\ \ldots \\ \ldots \\ {{Hp}\left( {{BrN}_{e}N_{r}} \right)} \end{pmatrix} = \begin{pmatrix} {{Hp}\left( {E1}\rightarrow{R1} \right)} \\ {{Hp}\left( {E1}\rightarrow{R2} \right)} \\ \ldots \\ \frac{{Hp}\left( {E1}\rightarrow R_{N_{r}} \right)}{{Hp}\left( {E2}\rightarrow{R1} \right)} \\ \ldots \\ \ldots \\ {{Hp}\left( E_{N_{e}}\rightarrow R_{N_{r}} \right)} \end{pmatrix}}$

The disrupted model Hp_(i) of the measurements is computed by a module 28 for computing the disrupted model based on the measurements Ip_(i) and the prediction H_(i) by the following equations:

${P_{i} = {\rho_{i}\frac{{jI}_{i}}{❘I_{i}❘}e^{j\alpha_{i}}}}{I_{i} = {H_{i} = {{Ip}_{i} - P_{i}}}}$

where P_(i) is the disruption made to the measurements Ip_(i) by the magnetic-field-disrupting material 3, ρ_(i) is the intensity of the disruption P_(i), α_(i) is a phase shift angle caused by the magnetic-field-disrupting material 3 in the disruption P_(i).

The phase shift angles α_(i) are related by the relationship characteristic of the magnetic-field-disrupting material 3 according to the equation α_(i)=−arctan(β·ω_(i)).

where β is a parameter characteristic of the magnetic-field-disrupting material 3.

The parameter β is computed as being identical for all the measurements Ip_(i) and represents a ratio of an inductance L_(p) modelling the disruptive material 3 to a resistance R_(p) modelling the disruptive material 3 in a certain position relative to the transmitter 10, the impedance of the disruptive material 3 being Zp=R_(p)+j·L_(p)·ω_(i).

The disruption P_(i) and/or the intensity ρ_(i) of the disruption P_(i) and/or the phase shift angle α_(i) and/or the parameter β and/or the disrupted model Hp_(i) is computed in such a way as to minimize the error E_(i).

The vector projection function Hp(H,Ip,β) is modeled to identify the parameter β.

For i ranging from 1 to N, N equations are to be solved.

Supposing that I_(i)=H_(i) and Ip_(i)=Hp_(i), this gives H _(i) +P _(i) −Ip _(i)=0

then the projected model vector Hp is obtained, the coordinates of which are Hp_(i) given by the following equation:

${{Hp}_{i} = {H_{i} \cdot \left( {1 + {\Delta_{i} \cdot \left( {j + {\beta \cdot \omega_{i}}} \right)}} \right)}}{with}{\Delta_{i} = \frac{\rho_{i}}{{❘H_{i}❘} \cdot \sqrt{{\beta^{2} \cdot \omega_{i}^{2}} + 1}}}$

Δ_(i) is called the projection distance.

The disruption P_(i) and/or the intensity ρ_(i) of the disruption P_(i) and/or the phase shift angle α_(i) and/or the parameter β and/or the disrupted model Hp_(i) and/or the projection distance Δ_(i) can be computed by the processing unit 25 by a computing method using successive iterations according to one of FIGS. 8 to 11. For example, according to an embodiment, at each iteration

-   -   the predictions H_(i) of the module 26 are initialized at given         values,     -   then the module 28 is used to compute by one of the described         computing methods or algorithms the disrupted model Hp_(i), as a         function of H_(i) and Ip_(i), for example according to the         equations         Δ_(i)=[         (H _(i))·         (Ip _(i))−         (H _(i))·         (Ip _(i))]/[(         (H _(i))²+         (H _(i))²]         β=[−         (H _(i))²−         (H _(i))²+         (H _(i))·         (Ip _(i))+         (H _(i))·         (Ip _(i))]/[(         (H _(i))·         (Ip _(i))−         (H _(i))·         (Ip _(i)))·ω_(i)],         Hp _(i) =H _(i)·(1+Δ_(i)·(j+β·ω _(i))),     -   then the module 27 is used to compute the error E_(i) according         to the equation         E _(i) =Ip _(i) −Hp _(i),     -   then the module 29 is used to compute the criterion C according         to the equation

${C = {\sum\limits_{i = 0}^{N}{❘{{HP}_{i} - {IP}_{i}}❘}^{2}}},$ or according to the equation

$C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}$

-   -   then the model of the module 26 is used to compute the         prediction H_(i) corresponding to the criterion C using the         minimization algorithm, for example of Levenberg-Marquardt type,         of the model. This minimization algorithm searches via the         module 26 for the position P and/or the orientation Q of the         receiver 20 which minimizes the deviation between the prediction         H_(i) according to the criterion C and the measurements Ip_(i).         The prediction H_(i) thus computed by the model 26 is used for         the following iteration. The iterations are continued until, for         example, the computed criterion C becomes less than a prescribed         non-zero positive bound η.

In the case where the disruptive material 3 is moving with respect to the transmitter 10 and/or the receiver 20, the parameter β varies over time as a function of the movement of the disruptive material 3.

Each module 26, 27, 28, 29 of the processing unit 25 and/or each step of the method can be implemented by computing means, for example by one or more calculator(s), one or more computers, one or more microprocessors and one or more computer programs. The computer program(s) comprise code instructions for implementing the compensating method.

According to an embodiment, the magnetic locator 1 comprising:

-   -   at least one transmitter 10 comprising a plurality (Ne) of         transmitter coils and a generator able to generate at least one         transmission signal, connected to the transmitter coils, so that         each transmitter coil is able to transmit at least one         transmitting magnetic field at at least one determined frequency         ω_(i) in response to the transmission signal that is sent to it         by the generator,     -   at least one receiver comprising one or more (Nr) receiver coils         and a measuring device, which is able to be connected to the         receiver coil or coils and which is able to measure for each         receiver coil a plurality of receiving magnetic fields         respectively induced by the plurality of transmitting magnetic         fields and called measurements Ip_(i) for i ranging from 1 to N,     -   a processing unit 25 comprising a field model making it possible         to compute a position P and/or an orientation Q of the receiver         by computing the prediction H_(i) of the measurements as a         function of the criterion C computed as a function of the error         E_(i), itself computed with respect to the measurements Ip_(i),

characterized in that

the processing unit 25 is configured so that the error E_(i) is computed by successive iterations from initial prescribed values of the prediction H_(i) as being the difference between the measurements Ip_(i) and a disrupted model Hp_(i) of the measurements, according to the equation E _(i) =Ip _(i) −Hp _(i),

the disrupted model Hp_(i) of the measurements verifying the following equations

${{Hp}_{i} = {H_{i} + P_{i}}},{P_{i} = {\rho_{i}\frac{{jH}_{i}}{❘H_{i}❘}e^{j\alpha_{i}}}},{\alpha_{i} = {- {\arctan\left( {\beta \cdot \omega_{i}} \right)}}},$

where P_(i) is a disruption made to the measurements Ip_(i) by the magnetic-field-disrupting material 3, ρ_(i) is the intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material 3, β being a parameter of the magnetic-field-disrupting material 3,

the parameter β being identical for all the measurements Ip_(i),

the processing unit 25 being configured to carry out the computation in such a way as to minimize the criterion C.

The processing unit 25 is thus configured to compensate for the presence of the magnetic-field-disrupting material 3.

The magnetic locator 1 thus comprises in its processing unit 25 means for compensating for the presence of the magnetic-field-disrupting material 3.

Below is a description of a first algorithm for carrying out the computation, with reference to the FIG. 9.

According to an embodiment, starting from the equation

H_(i)−Ip_(i)+Δ_(i)·H_(i)·(j+β·ω_(i)))=0, one obtains on the real parts and the imaginary parts thereof a system of two equations, which gives a unique solution for Δ_(i) and β, namely the solution Δ_(i)=[

(H _(i))·

(Ip _(i))−

(H _(i))·

(Ip _(i))]/[(

(H _(i))²+

(H _(i))²] and β=[−

(H _(i))²−

(H _(i))²+

(H _(i))·

(Ip _(i))+

(H _(i))·

(Ip _(i))]/[(

(H _(i))·

(Ip _(i))−

(H _(i))·

(Ip _(i)))·ω_(i)]

According to an embodiment, the first computing algorithm A makes provision for to be computed as being a solution to the following system of equations, resulting from the preceding equations giving Hp_(i), namely Δ_(i)·

(H _(i))·β·ω_(i)=

(Ip _(i))−

(H _(i))+Δ_(i)·

(H _(i)) and Δ_(i)·

(H _(i))·β·ω_(i)=

(Ip _(i))−

(H _(i))−Δ_(i)·

(H _(i)) for i ranging from 1 to N, where

(H_(i)) is the real part of the prediction H_(i) and is also written

(H_(i)),

(H_(i)) is the imaginary part of the prediction H_(i) and is also written ℑ(H_(i)),

(Ip_(i)) is the real part of the measurement Ip_(i) and is also written

(IP_(i))

(Ip_(i)) is the imaginary part of the measurement Ip_(i) and is also written ℑ(IP_(i)).

According to an embodiment, this first computing algorithm A can make provision for the parameter to be computed as minimizing a vector A·β−B over its coordinates, where

A is a vector, the first and second coordinates of which are respectively formed by: Δ_(i)·

(H_(i))·ω_(i) and Δ_(i)·

(H_(i))·ω_(i),

B is a vector, the first and second coordinates of which are respectively formed by:

(Ip_(i))−

(H_(i))+Δ_(i)·

(H_(i)) and

(Ip_(i))−

(H_(i))−Δ_(i)·

(H_(i)).

In an embodiment, the parameter β can be computed as minimizing the sum of the squares of the coordinates of the vector A·β−B for i ranging from 1 to N.

In another embodiment, the parameter β can be computed as being a minimal solution to the vector equation A·β=B, by

$\beta = \frac{A^{T} \cdot B}{{❘A❘}^{2}}$

The operator ^(T) denotes transposition.

In another embodiment, a weighting matrix denoted W is used such that: (W·A)·β=(W·B) with W=P·I and P=(p ₁ ,p ₂ ,p ₃ , . . . ,p _(N)).

The parameter β can then be computed as being

$\beta = \frac{\left( {W \cdot A} \right)^{T} \cdot B}{{❘\left( {W \cdot A} \right)❘}^{2}}$

For example, the p_(i) can be equal to a power of the norm of the H_(i), i.e. P=(|H ₀|^(k) ,|H ₁|^(k) , . . . ,|H _(N)|^(k)).

Thus, the greater k is, the more the large measurements are taken into account, in relation to the small ones. k=2 seems to be a good compromise. This avoids the small values of the coordinates of the vectors A and B leading to erroneous phases and in the estimating of β they can be rejected more often.

In an embodiment, this first computing algorithm A can be implemented by the following steps:

-   -   computing Δ_(i) according to the equation Δ_(i)=[         (H_(i))·         (Ip_(i))−         (H_(i))·         (Ip_(i))]/[(         (H_(i))²+         (H_(i))²],     -   computing β according to the equation

${\beta = {{\frac{A^{T} \cdot B}{{❘A❘}^{2}}{or}\beta} = \frac{\left( {W \cdot A} \right)^{T} \cdot B}{{❘\left( {W \cdot A} \right)❘}^{2}}}},$ computing Hp_(i) according to the equation Hp_(i)=H_(i)·(1+Δ_(i)·(j+β·ω_(i))), as represented in FIG. 9.

Below is a description of a second computing algorithm, with reference to FIG. 10.

According to an embodiment, starting from the equation H _(i) −Ip _(i)+Δ_(i) ·H _(i)·(j+β·ω _(i)))=0,

one obtains on the real parts and the imaginary parts thereof a system of two equations, which gives a unique solution for Δ_(i) and β, namely the solution Δ_(i)=[

(Ip _(i))−

(H _(i))]/[

(H _(i))+

(H _(i))·β·ω_(i)]

According to an embodiment, the second computing algorithm B makes provision for β to be computed as being a solution to the following equation, resulting from the preceding equation giving Δ_(i), namely β·ω_(i)·[

(H_(i))·

(Ip_(i))−

(Ip_(i))·

(H_(i))]+

(H_(i))²+

(H_(i))²−(

(H_(i))·

(Ip_(i))+

(H_(i))·(Ip_(i)))=0 for i ranging from 1 to N.

According to an embodiment, this second computing algorithm B can make provision for the parameter β to be computed as minimizing a vector A′·β+B′ over its coordinates, where

A′ is a vector, the coordinates of which are respectively formed by: ω_(i)·[

(H_(i))·

(Ip_(i))−

(Ip_(i))·

(H_(i))],

B′ is a vector, the coordinates of which are respectively formed by:

(H_(i))²+

(H_(i))²−(

(H_(i))·

(Ip_(i))+

(H_(i))·(Ip_(i))).

In an embodiment, the parameter β can be computed as minimizing the sum of the squares of the coordinates of the vector A′·β+B′ for i ranging from 1 to N.

In another embodiment, the parameter β can be computed as being a minimal solution to the vector equation A′·β=−B′, by

$\beta = \frac{A^{\prime T} \cdot B^{\prime}}{{❘A^{\prime}❘}^{2}}$

In an embodiment, this second computing algorithm B can be implemented by the following steps:

-   -   computing the vectors A′ and B′ according to the coordinates         indicated above,     -   computing β according to the equation

$\beta = \frac{A^{\prime T} \cdot B^{\prime}}{{❘A^{\prime}❘}^{2}}$ or according to the equation

$\beta = \frac{\left( {W \cdot A^{\prime}} \right)^{T} \cdot B^{\prime}}{{❘\left( {W \cdot A^{\prime}} \right)❘}^{2}}$ if the aforementioned weighting matrix W is used,

-   -   computing Δ_(i) according to the equation Δ_(i)=[         (Ip_(i))−         (H_(i))]/[         (H_(i))+         (H_(i))·β·ω_(i)] and computing Hp_(i) according to the equation         Hp_(i)=H_(i)·(1+Δ_(i)·(j+P·ω_(i))), as represented in FIG. 10.

Below is a description of a third computing algorithm C, with reference to FIG. 11.

According to an embodiment, starting from the equation H _(i) −Ip _(i)+Δ_(i) ·H _(i)·(j+β·ω _(i)))=0,

one obtains on the real parts and the imaginary parts thereof a system of two equations namely

(H _(i))−(Ip _(i))+Δ_(i)·(

(H _(i))·β·ω_(i)−

(H _(i)))=0 and

(H _(i))−

(Ip _(i))+Δ_(i)·(

(H _(i))+

(H _(i))·β·ω_(i))=0 According to an embodiment, the third computing algorithm C makes provision for β to be computed as being a solution to this system of equations. For example the Levenberg-Marquardt algorithm, which produces the inverse solution, also computes the parameter β as being a state variable.

According to an embodiment, this third computing algorithm C can make provision for a vector Δ having Δ_(i) as coordinates to be computed as minimizing a vector C·Δ−D over its coordinates, where

D is a vector, the first and second coordinates of which are respectively formed by: −(

(H_(i))−

(Ip_(i))) and −(

(H_(i))−

(Ip_(i))),

C is a matrix, having as coefficients corresponding to the Δ_(i) respectively

(H_(i))·β·ω_(i)−

(H_(i)) and

(H_(i))+

(H_(i))·β·ω_(i), i.e. the vector C·Δ−D is equal to:

${\begin{pmatrix} \left( {{\left( H_{0} \right) \cdot \beta \cdot \omega_{0}} - \left( H_{0} \right)} \right) & 0 & 0 \\ \left( {\left( H_{0} \right) + {\left( H_{0} \right) \cdot \beta \cdot \omega_{0}}} \right) & 0 & 0 \\ 0 & \left( {{\left( H_{1} \right) \cdot \beta \cdot \omega_{1}} - \left( H_{1} \right)} \right) & 0 \\ 0 & \left( {\left( H_{1} \right) + {\left( H_{1} \right) \cdot \beta \cdot \omega_{1}}} \right) & 0 \\ \ldots & \ldots & \ldots \\ 0 & 0 & \left( {{\left( H_{N} \right) \cdot \beta \cdot \omega_{N}} - \left( H_{N} \right)} \right) \\ 0 & 0 & \left( {\left( H_{N} \right) + {\left( H_{N} \right) \cdot \beta \cdot \omega_{N}}} \right) \end{pmatrix}\begin{pmatrix} \Delta_{0} \\ \Delta_{1} \\ \ldots \\ \Delta_{N} \end{pmatrix}} - \begin{pmatrix} {\left( {IP}_{0} \right) - \left( H_{0} \right)} \\ {\left( {IP}_{0} \right) - \left( H_{0} \right)} \\ {\left( {IP}_{1} \right) - \left( H_{1} \right)} \\ {\left( {IP}_{1} \right) - \left( H_{1} \right)} \\ \ldots \\ \ldots \\ {\left( {IP}_{N} \right) - \left( H_{N} \right)} \\ {\left( {IP}_{N} \right) - \left( H_{N} \right)} \end{pmatrix}$

In another embodiment, the parameter β can be computed as being a solution to the vector equation C·Δ=D. For example, a matrix C⁺ is computed which is the pseudoinverse of C, to compute the vector Δ as being Δ=C⁺·D.

In an embodiment, this third computing algorithm C can be implemented by the following steps:

-   -   estimating the parameter β via the Levenberg-Marquardt         algorithm,     -   computing the matrix C and the vector D, computing the vector Δ         verifying the equation C·Δ=D, for example by Δ=C⁺·D     -   computing Hp_(i) according to the equation         Hp_(i)=H_(i)·(1+Δ_(i)·(j+P·ω_(i))), as represented in FIG. 11.

Below is a description of a variant of the compensating method and device according to the invention. According to this variant, each transmitter coil transmits several frequencies. Thus several frequencies are transmitted per transmitter coil axis. It can thus be considered that the magnetic-field-disrupting material 3 constitutes several magnetic-field-disrupting materials 3 and that the projection distances are linked for the different frequencies. It is for example possible to consider a locator 1 having a transmitter coil transmitting at K different frequencies and J receiver coils, with K greater than or equal to 2 and J greater than or equal to 2. This gives N=K·J measurements Ip_(i). This gives N equations of the form: Ip _(i) =H _(i)·(1+Δ_(i)·(j+β·ω _(i))).

It is possible to break down the index i belonging to [1, N] into two indices j belonging to [1, J] and k belonging to [1, K]. The preceding equation then becomes: Ip _(jk) =H _(jk)·(1+Δ_(jk)·(j+β·ω _(jk))).

In all the equations, the letter j not appearing in the index is the complex number, the square of which is equal to −1, as is known to those skilled in the art.

The change of variable is set:

Δ jk = ρ j ⁢ ω j ⁢ k 1 + ( β · ω jk ) 2

ρ_(j) depends only on j.

The equations to be minimized become:

$F_{({\rho_{j},\beta})} = \left\{ \begin{matrix} \begin{matrix} {R_{jk} = {\left( H_{jk} \right) - \left( {Ip}_{jk} \right) + {\rho_{j}{\frac{\omega_{jk}}{1 + \left( {\beta{\cdot \omega_{jk}}} \right)^{2}} \cdot}}}} \\ {\left( {{\left( H_{jk} \right) \cdot \beta \cdot \omega_{jk}} - \left( H_{jk} \right)} \right) = 0} \end{matrix} \\ \begin{matrix} {I_{jk} = {\left( H_{jk} \right) - \left( {Ip}_{jk} \right) + {\rho_{j}{\frac{\omega_{jk}}{1 + \left( {\beta{\cdot \omega_{jk}}} \right)^{2}} \cdot}}}} \\ {\left( {{\left( H_{jk} \right) \cdot \beta \cdot \omega_{jk}} + \left( H_{jk} \right)} \right) = 0} \end{matrix} \end{matrix} \right.$

This gives a system of 2JK equations with J+1 unknowns to be minimized.

According to an embodiment, 2·J·K≥J+1.

Solving two by two obtains: v _(jk)=

(H _(jk))·

(Ip _(jk))+ℑ(H _(jk))·ℑ(Ip _(jk)) u _(jk)=(

(H _(jk))·ℑ(Ip _(jk))−ℑ(H _(jk))·

(Ip _(jk)))·ω_(jk) u _(jk)·ρ_(j) =|H _(jk)|² +|Ip _(jk)|²−2·v _(jk) u _(jk) ·β=v _(jk) −|H _(jk)|²

The least squares solution is expressed by:

${\rho_{j} = \frac{\sum\limits_{k = 1}^{K}{u_{jk} \cdot \left( {{❘H_{jk}❘}^{2} + {❘{Ip}_{jk}❘}^{2} - {2 \cdot v_{jk}}} \right)}}{\sum\limits_{k = 1}^{K}u_{jk}^{2}}}{and}{\beta = \frac{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}{u_{jk} \cdot \left( {v_{jk} - {❘H_{jk}❘}^{2}} \right)}}}{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}u_{jk}^{2}}}}$

If there are L transmitters 10 with L greater than or equal to 2, there will be a third index l belonging to [l, L] and one parameter β_(l) is searched for per transmitter 10 _(l) according to the following equations:

${\rho_{jl} = \frac{\sum\limits_{k = 1}^{K}{u_{jkl} \cdot \left( {{❘H_{jkl}❘}^{2} + {❘{Ip}_{jkl}❘}^{2} - {2 \cdot v_{jkl}}} \right)}}{\sum\limits_{k = 1}^{K}u_{jkl}^{2}}}{and}{\beta_{l} = \frac{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}{u_{jkl} \cdot \left( {v_{jkl} - {❘H_{jkl}❘}^{2}} \right)}}}{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}u_{jkl}^{2}}}}$

In this case there are 2JKL measurements for L(J+1) unknowns. According to an embodiment, 2·J·K·L≥L·(J+1).

It is found that increasing the number of frequencies per axis does not increase the number of unknowns and therefore increases the redundancy.

Of course, the possibilities, embodiments, features, variants and examples above can be combined with one another or be selected independently of one another. 

The invention claimed is:
 1. A method for compensating a magnetic locator in the presence of at least one magnetic-field-disrupting material, said method performed by a magnetic locator comprising: at least one transmitter comprising at least one transmitter coil and at least one generator of at least one transmission signal, connected to the at least one transmitter coil, so that the at least one transmitter coil transmits at least one transmitting magnetic field at at least one determined angular frequency ω_(i) in response to the transmission signal that is sent to it by the generator, at least one receiver comprising at least one receiver coil and a measuring device, which is connected to the at least one receiver coil and which supplies at least one measurement Ip_(i) of a receiving magnetic field respectively induced by the transmitting magnetic field in each receiver coil, in such a way as to supply the at least one measurement Ip_(i) for i ranging from 1 to N, with N=Nr·Ne, wherein Ne is a first number of transmitter coil, Nr is a second number of receiver coil, Ne≥1, Nr≥1, a processing unit, wherein the method comprises the followings steps: initializing by a first calculator of the processing unit a prediction H_(i) of the at least one measurement at prescribed initial values, computing in successive iterations by a second calculator of the processing unit a criterion C until the criterion C becomes less than a prescribed non-zero positive bound η, wherein at each iteration: computing by a third calculator of the processing unit a value Δ_(i), a parameter β of the magnetic-field-disrupting material so as to minimize a vector A·β−B over its coordinates, and a disrupted model Hp_(i) as a function of H_(i) and Ip_(i) according to the equations Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢e j ⁢ α i , α i = - arctan ⁡ ( β · ω i ) , where P_(i) is a disruption made to the at least one measurement Ip_(i) by the magnetic-field-disrupting material, ρ_(i) is an intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, the parameter β being identical for all the measurements Ip_(i), Δi=[R(H _(i))·I(Ip _(i))−I(H _(i))·R(Ip _(i))]/[I(H _(i))² +R(H _(i))²], Hp _(i) =H _(i)·(1+Δ_(i)·(j+β·ω _(i))), where A is a vector, first and second coordinates of which are respectively formed by: A_(2i)=Δ_(i)·R(H_(i))·ω_(i) and A_(2i+1)=Δ_(i)·I(H_(i))·ω_(i), B is a vector, first and second coordinates of which are respectively formed by: B _(2i) =R(Ip _(i))−R(H _(i))+Δ_(i) ·I(H _(i)) and B _(2i+1) =I(Ip _(i))−I(H _(i))−Δ_(i) ·R(H _(i)), for i ranging from 1 to N, where R(H_(i)) is a real part of the prediction H_(i), I(H_(i)) is an imaginary part of the prediction H_(i), R(Ip_(i)) is a real part of the at least one measurement Ip_(i), I(Ip_(i)) is an imaginary part of the at least one measurement Ip_(i), then computing by a fourth calculator of the processing unit the error E according to the equation E _(i) =Ip _(i) −Hp _(i), computing by the second calculator the criterion C according to the equation ${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$ then computing by a field model of the first calculator at least one of a position of the at least one receiver and of an orientation of the at least one receiver by computing the prediction H_(i) corresponding to the criterion C.
 2. A method for compensating a magnetic locator in the presence of at least one magnetic-field-disrupting material, said method performed by a magnetic locator comprising: at least one transmitter comprising at least one transmitter coil and at least one generator of at least one transmission signal, connected to the at least one transmitter coil, so that the at least one transmitter coil transmits at least one transmitting magnetic field at at least one determined angular frequency ω_(i) in response to the transmission signal that is sent to it by the generator, at least one receiver comprising at least one receiver coil and a measuring device, which is connected to the at least one receiver coil and which supplies at least one measurement Ip_(i) of a receiving magnetic field respectively induced by the transmitting magnetic field in each receiver coil, in such a way as to supply the at least one measurement Ip_(i) for i ranging from 1 to N, with N=Nr·Ne, wherein Ne is a first number of transmitter coil, Nr is a second number of receiver coil, Ne≥1, Nr≥1, a processing unit, wherein the method comprises the followings steps: initializing by a first calculator of the processing unit a prediction H_(i) of the at least one measurement at prescribed initial values, computing in successive iterations by a second calculator of the processing unit a criterion C until the criterion C becomes less than a prescribed non-zero positive bound η, wherein at each iteration: computing by a third calculator of the processing unit a value Δ_(i), a parameter β of the magnetic-field-disrupting material and a disrupted model Hp_(i) as a function of H_(i) and Ip_(i) according to the equations Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢ e j ⁢ α i , α i = - arctan ⁡ ( β · ω i ) , where P_(i) is a disruption made to the at least one measurement Ip_(i) by the magnetic-field-disrupting material, ρ_(i) is an intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, β being a parameter of the magnetic-field-disrupting material, the parameter β being identical for all the measurements Ip_(i), ${{\Delta_{i} = {\left\lbrack {{R\left( H_{i} \right)}{{\cdot {I\left( {Ip}_{i} \right)}} - {{I\left( H_{i} \right)} \cdot {R\left( {Ip}_{i} \right)}}}} \right\rbrack/\left\lbrack \left( {{I\left( H_{i} \right)}^{2} + {R\left( H_{i} \right)}^{2}} \right. \right\rbrack}},{\beta = \frac{A^{T} \cdot B}{{❘A❘}^{2}}}}{{{Hp}_{i} = {H_{i} \cdot \left( {1 + {\Delta_{i} \cdot \left( {j + {\beta \cdot \omega_{i}}} \right)}} \right)}},}$ where A is a vector, first and second coordinates of which are respectively formed by: A_(2i)=Δ_(i)·R(H_(i))·ω_(i) and A_(2i+1)=Δ_(i)·I(H_(i))·ω_(i), B is a vector, first and second coordinates of which are respectively formed by: B_(2i)=R(Ip_(i))−R(H_(i))+Δ_(i)·I(H_(i)) and B_(2i+1)=I(Ip_(i))−I(H_(i))−Δ_(i)·R(H_(i)), for i ranging from 1 to N, where R(H_(i)) is a real part of the prediction H_(i), I(H_(i)) is an imaginary part of the prediction H_(i), R(Ip_(i)) is a real part of the at least one measurement Ip_(i), I(Ip_(i)) is an imaginary part of the at least one measurement Ip_(i), then computing by a fourth calculator of the processing unit the error E_(i) according to the equation E _(i) =Ip _(i) −Hp _(i), computing by the second calculator the criterion C according to the equation ${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$ then computing by a field model of the first calculator at least one of a position of the at least one receiver and of an orientation of the at least one receiver by computing the prediction H_(i) corresponding to the criterion C.
 3. A method for compensating a magnetic locator in the presence of at least one magnetic-field-disrupting material, said method performed by a magnetic locator comprising: at least one transmitter comprising at least one transmitter coil and at least one generator of at least one transmission signal, connected to the at least one transmitter coil, so that the at least one transmitter coil transmits at least one transmitting magnetic field at at least one determined angular frequency ω_(i) in response to the transmission signal that is sent to it by the generator, at least one receiver comprising at least one receiver coil and a measuring device, which is connected to the at least one receiver coil and which supplies at least one measurement Ip_(i) of a receiving magnetic field respectively induced by the transmitting magnetic field in each receiver coil, in such a way as to supply the at least one measurement Ip_(i) for i ranging from 1 to N, with N=Nr·Ne, wherein Ne is a first number of transmitter coil, Nr is a second number of receiver coil, Ne≥1, Nr≥1, a processing unit, wherein the method comprises the followings steps: initializing by a first calculator of the processing unit a prediction H_(i) of the at least one measurement at prescribed initial values, computing in successive iterations by a second calculator of the processing unit a criterion C until the criterion C becomes less than a prescribed non-zero positive bound η, wherein at each iteration: computing by a third calculator of the processing unit a parameter β of the magnetic-field-disrupting material so as to minimize a vector A′·β+B′ over its coordinates, then Δ_(i) and a disrupted model Hp_(i) as a function of H_(i) and Ip_(i) according to the equations: Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢ e j ⁢ α i , α i = - arctan ⁡ ( β · ω i ) , where P_(i) is a disruption made to the at least one measurement Ip_(i) by the magnetic-field-disrupting material, ρ_(i) is an intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, the parameter β being identical for all the measurements Ip_(i), Δ_(i)=(I(Ip _(i))−I(H _(i)))/(R(H _(i))+I(H _(i))βω_(i)) Hp _(i) =H _(i)·(1+Δ_(i)·(j+β·ω _(i))) where A′ is a vector, coordinates of which are respectively formed by: ω_(i)·[R(H _(i))·I(Ip _(i))−R(Ip _(i))·I(H _(i))], B′ is a vector, the coordinates of which are respectively formed by: R(H _(i))² +I(H _(i))²−(R(H _(i))·R(Ip _(i))+I(H _(i))·(Ip _(i))), for i ranging from 1 to N, where R(H_(i)) is a real part of the prediction H_(i), I(H_(i)) is an imaginary part of the prediction H_(i), R(Ip_(i)) is a real part of the at least one measurement Ip_(i), I(Ip_(i)) is an imaginary part of the at least one measurement Ip_(i), then computing by a fourth calculator of the processing unit the error E_(i) according to the equation E _(i) =Ip _(i) −Hp _(i), computing by the second calculator the criterion C according to the equation ${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$ then computing by a field model of the first calculator at least one of a position of the at least one receiver and of an orientation of the at least one receiver by computing the prediction H_(i) corresponding to the criterion C.
 4. A method for compensating a magnetic locator in the presence of at least one magnetic-field-disrupting material, said method performed by a magnetic locator comprising: at least one transmitter comprising at least one transmitter coil and at least one generator of at least one transmission signal, connected to the at least one transmitter coil, so that the at least one transmitter coil transmits at least one transmitting magnetic field at at least one determined angular frequency ω_(i) in response to the transmission signal that is sent to it by the generator, at least one receiver comprising at least one receiver coil and a measuring device, which is connected to the at least one receiver coil and which supplies at least one measurement Ip_(i) of a receiving magnetic field respectively induced by the transmitting magnetic field in each receiver coil, in such a way as to supply the at least one 1 measurement Ip_(i) for i ranging from 1 to N, with N=Nr·Ne, wherein Ne is a first number of transmitter coil, Nr is a second number of receiver coil, Ne≥1, Nr≥1, a processing unit, wherein the method comprises the followings steps: initializing by a first calculator of the processing unit a prediction H_(i) of the at least one measurement at prescribed initial values, computing in successive iterations by a second calculator of the processing unit a criterion C until the criterion C becomes less than a prescribed non-zero positive bound η, wherein at each iteration: computing by a third calculator of the processing unit a parameter β then a value Δ_(i) and a disrupted model Hp_(i) as a function of H_(i) and Ip_(i) according to the equations: Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢ e j ⁢ α i , α i = - arctan ⁡ ( β · ω i ) , where P_(i) is a disruption made to the at least one measurement Ip_(i) by the magnetic-field-disrupting material, ρ_(i) is an intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, the parameter β being identical for all the measurements Ip_(i), ${\beta = \frac{A^{\prime T} \cdot B^{\prime}}{{❘A^{\prime}❘}^{2}}}{{\Delta_{i =}\left( {{I\left( {Ip}_{i} \right)} - {I\left( H_{i} \right)}} \right)}/\left( {{R\left( H_{i} \right)} + {{I\left( H_{i} \right)}\beta\omega_{i}}} \right)}{{Hp}_{i} = {H_{i} \cdot \left( {1 + {\Delta_{i} \cdot \left( {j + {\beta \cdot \omega_{i}}} \right)}} \right)}}$ where A′ is a vector, coordinates of which are respectively formed by: ω_(i)·[R(H _(i))·I(Ip _(i))−R(Ip _(i))·I(H _(i))], B′ is a vector, coordinates of which are respectively formed by: R(H _(i))² +I(H _(i))²−(R(H _(i))·R(Ip _(i))+I(H _(i))·(Ip _(i))), for i ranging from 1 to N, where R(H_(i)) is a real part of the prediction H_(i), I(H_(i)) is an imaginary part of the prediction H_(i), R(Ip_(i)) is a real part of the at least one measurement Ip_(i), I(Ip_(i)) is an imaginary part of the at least one measurement Ip_(i), then computing by a fourth calculator of the processing unit the error E_(i) according to the equation E _(i) =Ip _(i) −Hp _(i), computing by the second calculator the criterion C according to the equation ${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$ then computing by a field model of the first calculator at least one of a position of the at least one receiver and of an orientation of the at least one receiver by computing the prediction H_(i) corresponding to the criterion C.
 5. The method as claimed in any one of claims 1 to 4, wherein the field model to compute the position and/or the orientation of the at least one receiver by computing the prediction H_(i) of the at least one measurement as a function of the criterion C uses a Levenberg-Marquardt minimization algorithm.
 6. The method as claimed in any one of claims 1 to 4, wherein the determined angular frequencies ω_(i) are separate from one another.
 7. The method as claimed in any one of claims 1 to 4, wherein the at least one transmitter comprises a plurality of transmitter coils respectively transmitting a plurality of transmitting magnetic fields, with Ne≥2 the measuring device measuring for each receiver coil a plurality of receiving magnetic fields respectively induced by the plurality of transmitting magnetic fields in the receiver coil and forming the measurements Ip_(i) for i ranging from 1 to N, with N≥2.
 8. The method as claimed in any one of the claims 1 to 4, wherein the magnetic locator comprises a transmitter coil transmitting at K different angular frequencies ω_(jk) and J receiver coils, the measuring device supplying the measurements Ip_(jk) for the index j ranging from 1 to J and the index k ranging from 1 to K, the processing unit computes ${\rho_{j} = \frac{\sum\limits_{k = 1}^{K}{u_{jk} \cdot \left( {{❘H_{jk}❘}^{2} + {❘{Ip}_{jk}❘}^{2} - {2 \cdot v_{jk}}} \right)}}{\sum\limits_{k = 1}^{K}u_{jk}^{2}}}{and}{\beta = \frac{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}{u_{jk} \cdot \left( {v_{jk} - {❘H_{jk}❘}^{2}} \right)}}}{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}u_{jk}^{2}}}}{with}{v_{jk} = {{{R\left( H_{jk} \right)} \cdot {R\left( {Ip}_{jk} \right)}} + {{I\left( H_{jk} \right)} \cdot {I\left( {Ip}_{jk} \right)}}}}{u_{jk} = {\left( {{{R\left( H_{jk} \right)} \cdot {I\left( {Ip}_{jk} \right)}} - {{I\left( H_{jk} \right)} \cdot {R\left( {Ip}_{jk} \right)}}} \right) \cdot \omega_{jk}}}{{Ip}_{jk} = {H_{jk} \cdot \left( {1 + {\Delta_{jk} \cdot \left( {j + {\beta \cdot \omega_{jk}}} \right)}} \right)}}{\Delta_{jk} = {\rho_{j}\frac{\omega_{jk}}{1 + \left( {\beta \cdot \omega_{jk}} \right)^{2}}}}$ where H_(jk) is the prediction, ρ_(j) is the intensity of the disruption, Hp_(jk) is the disrupted model and Hp_(jk)−Ip_(jk)=0.
 9. The method as claimed in any one of the claims 1-4, wherein the magnetic locator comprises L transmitter coils or L transmitters, which transmit at K different angular frequencies ω_(jkl), and J receiver coils, the measuring device supplying the measurements Ip_(jkl) for the index j ranging from 1 to J, the index k ranging from 1 to K and the index I ranging from 1 to L, the processing unit computes ${\rho_{jl} = \frac{\sum\limits_{k = 1}^{K}{u_{jkl} \cdot \left( {{❘H_{jkl}❘}^{2} + {❘{Ip}_{jkl}❘}^{2} - {2 \cdot v_{jkl}}} \right)}}{\sum\limits_{k = 1}^{K}u_{jkl}^{2}}}{and}{\beta_{l} = \frac{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}{u_{jkl} \cdot \left( {v_{jkl} - {❘H_{jkl}❘}^{2}} \right)}}}{\sum\limits_{j = 1}^{J}{\sum\limits_{k = 1}^{K}u_{jkl}^{2}}}}{with}{v_{jkl} = {{{R\left( H_{jkl} \right)} \cdot {R\left( {Ip}_{jkl} \right)}} + {{I\left( H_{jkl} \right)}{I\left( {Ip}_{jkl} \right)}}}}{u_{jkl} = {\left( {{{R\left( H_{jkl} \right)} \cdot {I\left( {Ip}_{jkl} \right)}} - {{I\left( H_{jkl} \right)} \cdot {R\left( {Ip}_{jkl} \right)}}} \right) \cdot \omega_{jkl}}}{{Ip}_{jkl} = {H_{jkl} \cdot \left( {1 + {\Delta_{jkl} \cdot \left( {j + {\beta \cdot \omega_{jkl}}} \right)}} \right)}}{\Delta_{jkl} = {\rho_{j}\frac{\omega_{jkl}}{1 + \left( {\beta \cdot \omega_{jkl}} \right)^{2}}}}$ where H_(jkl) is the prediction, ρ_(jl) is the intensity of the disruption, Hp_(jkl) is the disrupted model and Hp_(jkl)−Ip_(jkl)=0.
 10. A non-transitory computer-readable storage medium, comprising code instructions stored thereon that, when executed by a calculator, cause the calculator to execute the method for compensating a magnetic locator in the presence of at least one magnetic-field-disrupting material as claimed in any one of claims 1 to
 4. 11. A magnetic locator comprising: at least one transmitter comprising at least one transmitter coil and at least one generator able to generate at least one transmission signal, connected to the at least one transmitter coil, so that the at least one transmitter coil is able to transmit at least one transmitting magnetic field at at least one determined angular frequency ω_(i) in response to the transmission signal that is sent to it by the generator, at least one receiver comprising at least one receiver coil and a measuring device, which is able to be connected to the at least one receiver coil and which is able to supply at least one measurement Ip_(i) of a receiving magnetic field respectively induced by the transmitting magnetic field in each receiver coil, to supply the at least one measurement Ip_(i) for i ranging from 1 to N, with N=Nr·Ne, wherein Ne is a first number of transmitter coil, Nr is a second number of receiver coil, Ne≥1, Nr≥1, a processing unit comprising a first calculator for initializing a prediction H_(i) of the at least one measurement at prescribed initial values, a second calculator for computing in successive iterations a criterion C according to the equation ${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$ until the criterion C becomes less than a prescribed non-zero positive bound η, a third calculator for computing at each iteration a value Δ_(i), a parameter β of the magnetic-field-disrupting material so as to minimize a vector A·β−B over its coordinates, and a disrupted model Hp_(i) as a function of H_(i) and Ip_(i) according to the equations: Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢ e j ⁢ α i ⁢ α i = - arctan ⁡ ( β · ω i ) , where P_(i) is a disruption made to the at least one measurement Ip_(i) by the magnetic-field-disrupting material, ρ_(i) is an intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, the parameter β being identical for all the measurements Ip_(i), Δ_(i) = [R(H_(i)) ⋅ I(Ip_(i)) − I(H_(i)) ⋅ R(Ip_(i))]/[(I(H_(i))² + R(H_(i))²], Hp_(i) = H_(i) ⋅ (1 + Δ_(i) ⋅ (j + β ⋅ ω_(i))), where A is a vector, first and second coordinates of which are respectively formed by: A _(2i)=Δ_(i) ·R(H _(i))·ω_(i) and A _(2i+1)=Δ_(i) ·I(H _(i))·ω_(i), B is a vector, first and second coordinates of which are respectively formed by: B _(2i) =R(Ip _(i))−R(H _(i))+Δ_(i) ·I(H _(i)) and B _(2i+1) =I(Ip _(i))−I(H _(i))−Δ_(i) ·R(H _(i)), for i ranging from 1 to N, where R(H_(i)) is a real part of the prediction H_(i), I(H_(i)) is an imaginary part of the prediction H_(i), R(Ip_(i)) is a real part of the at least one measurement Ip_(i), I(Ip_(i)) is an imaginary part of the at least one measurement Ip_(i), a fourth calculator for computing at each iteration the error E_(i) according to the equation E _(i) =Ip _(i) −Hp _(i), the first calculator having a field model for computing at each iteration at least one of a position of the at least one receiver and of an orientation of the at least one receiver by computing the prediction H_(i) corresponding to the criterion C.
 12. A magnetic locator comprising: at least one transmitter comprising at least one transmitter coil and at least one generator able to generate at least one transmission signal, connected to the at least one transmitter coil, so that the at least one transmitter coil is able to transmit at least one transmitting magnetic field at at least one determined angular frequency ω_(i) in response to the transmission signal that is sent to it by the generator, at least one receiver comprising at least one receiver coil and a measuring device, which is able to be connected to the at least one receiver coil and which is able to supply at least one measurement Ip_(i) of a receiving magnetic field respectively induced by the transmitting magnetic field in each receiver coil, to supply the at least one measurement Ip_(i) for i ranging from 1 to N, with N=Nr·Ne, wherein Ne is a first number of transmitter coil, Nr is a second number of receiver coil, Ne≥1, Nr≥1, a processing unit comprising a first calculator for initializing a prediction H_(i) of the at least one measurement at prescribed initial values, a second calculator for computing in successive iterations a criterion C according to the equation ${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$ until the criterion C becomes less than a prescribed non-zero positive bound η, a third calculator for computing at each iteration a value Δ_(i), a parameter β of the magnetic-field-disrupting material and a disrupted model Hp_(i) as a function of H_(i) and Ip_(i) according to the equations Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢ e j ⁢ α i ⁢ α i = - arctan ⁡ ( β · ω i ) , where P_(i) is a disruption made to the at least one measurement Ip_(i) by the magnetic-field-disrupting material, ρ_(i) is an intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, the parameter β being identical for all the measurements Ip_(i), $\Delta_{i} = {\left\lbrack {{{R\left( H_{i} \right)} \cdot {I\left( {Ip}_{i} \right)}} - {{I\left( H_{i} \right)} \cdot {R\left( {Ip}_{i} \right)}}} \right\rbrack/\left\lbrack {\left( {{I\left( H_{i} \right)}^{2} + {R\left( H_{i} \right)}^{2}} \right\rbrack,{\beta = {\frac{A^{T} \cdot B}{{A}^{2}}{{Hp}_{i} = {H_{i} \cdot \left( {1 + {\Delta_{i} \cdot \left( {j + {\beta \cdot \omega_{i}}} \right)}} \right)}}}},} \right.}$ where A is a vector, first and second coordinates of which are respectively formed by: A_(2i)=Δ_(i)·R(H_(i))·ω_(i) and A_(2i+1)=Δ_(i)·I(H_(i))·ω_(i), B is a vector, first and second coordinates of which are respectively formed by: B _(2i) =R(Ip _(i))−R(H _(i))+Δ_(i) ·I(H _(i)) and B _(2i+1) =I(Ip _(i))−I(H _(i))−Δ_(i) ·R(H _(i)), for i ranging from 1 to N, where R(H_(i)) is a real part of the prediction H_(i), I(H_(i)) is an imaginary part of the prediction H_(i), R(Ip_(i)) is a real part of the at least one measurement Ip_(i), I(Ip_(i)) is an imaginary part of the at least one measurement Ip_(i), a fourth calculator for computing at each iteration the error E_(i) according to the equation E _(i) =Ip _(i) −Hp _(i), the first calculator having a field model for computing at each iteration at least one of a position of the at least one receiver and of an orientation of the at least one receiver by computing the prediction H_(i) corresponding to the criterion C.
 13. A magnetic locator comprising: at least one transmitter comprising at least one transmitter coil and at least one generator able to generate at least one transmission signal, connected to the at least one transmitter coil, so that the at least one transmitter coil is able to transmit at least one transmitting magnetic field at at least one determined angular frequency ω_(i) in response to the transmission signal that is sent to it by the generator, at least one receiver comprising at least one receiver coil and a measuring device, which is able to be connected to the at least one receiver coil and which is able to supply at least one measurement Ip_(i) of a receiving magnetic field respectively induced by the transmitting magnetic field in each receiver coil, to supply the at least one measurement Ip_(i) for i ranging from 1 to N, with N=Nr·Ne, wherein Ne is a first number of transmitter coil, Nr is a second number of receiver coil, Ne≥1, Nr≥1, a processing unit comprising a first calculator for initializing a prediction H_(i) of the at least one measurement at prescribed initial values, a second calculator for computing in successive iterations a criterion C according to the equation ${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$ until the criterion C becomes less than a prescribed non-zero positive bound η, a third calculator for computing at each iteration a parameter β of the magnetic-field-disrupting material so as to minimize a vector A′·β+B′ over its coordinates, then a value Δ_(i) and a disrupted model Hp_(i) as a function of H_(i) and Ip_(i) according to the equations: Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢ e j ⁢ α i ⁢ α i = - arctan ⁡ ( β · ω i ) , where P_(i) is a disruption made to the at least one measurement Ip_(i) by the magnetic-field-disrupting material, ρ_(i) is an intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, the parameter β being identical for all the measurements Ip_(i), Δ_(i)=(I(Ip _(i))−I(H _(i)))/(R(H _(i))+I(H _(i))βω_(i)) Hp _(i) =H _(i)·(1+Δ_(i)·(j+β·ω _(i))) where A′ is a vector, coordinates of which are respectively formed by: ω_(i)·[R(H _(i))·I(Ip _(i))−R(Ip _(i))·I(H _(i))], B′ is a vector, the coordinates of which are respectively formed by: R(H _(i))² +I(H _(i))²−(R(H _(i))·R(Ip _(i))+I(H _(i))·(Ip _(i))), for i ranging from 1 to N, where R(H_(i)) is a real part of the prediction H_(i), I(H_(i)) is an imaginary part of the prediction H_(i), R(Ip_(i)) is a real part of the at least one measurement Ip_(i), I(Ip_(i)) is an imaginary part of the at least one measurement Ip_(i), a fourth calculator for computing at each iteration the error E_(i) according to the equation E _(i) =Ip _(i) −Hp _(i), the first calculator having a field model for computing at each iteration at least one of a position of the at least one receiver and of an orientation of the at least one receiver by computing the prediction H_(i) corresponding to the criterion C.
 14. A magnetic locator comprising: at least one transmitter comprising at least one transmitter coil and at least one generator able to generate at least one transmission signal, connected to the at least one transmitter coil, so that the at least one transmitter coil is able to transmit at least one transmitting magnetic field at at least one determined angular frequency ω_(i) in response to the transmission signal that is sent to it by the generator, at least one receiver comprising at least one receiver coil and a measuring device, which is able to be connected to the at least one receiver coil and which is able to supply at least one measurement Ip_(i) of a receiving magnetic field respectively induced by the transmitting magnetic field in each receiver coil, to supply the at least one measurement Ip_(i) for i ranging from 1 to N, with N=Nr·Ne, wherein Ne is a first number of transmitter coil, Nr is a second number of receiver coil, Ne≥1, Nr≥1, a processing unit comprising a first calculator for initializing a prediction H_(i) of the at least one measurement at prescribed initial values, a second calculator for computing in successive iterations a criterion C according to the equation ${C = {\sum\limits_{i = 0}^{N}{❘E_{i}❘}^{2}}},$ until the criterion C becomes less than a prescribed non-zero positive bound η, a third calculator for computing at each iteration a parameter β of the magnetic-field-disrupting material then a value Δ_(i) and a disrupted model Hp_(i) as a function of H_(i) and Ip_(i) according to the equations: Hp i = H i + P i , P i = ρ i ⁢ jH i ❘ "\[LeftBracketingBar]" H i ❘ "\[RightBracketingBar]" ⁢ e j ⁢ α i ⁢ α i = - arctan ⁡ ( β · ω i ) , where P_(i) is a disruption made to the at least one measurement Ip_(i) by the magnetic-field-disrupting material, ρ_(i) is an intensity of the disruption, α_(i) is a phase shift angle caused by the magnetic-field-disrupting material, the parameter β being identical for all the measurements Ip_(i), ${\beta = \frac{A^{\prime T} \cdot B^{\prime}}{{❘A^{\prime}❘}^{2}}}{{\Delta_{i =}\left( {{I\left( {Ip}_{i} \right)} - {I\left( H_{i} \right)}} \right)}/\left( {{R\left( H_{i} \right)} + {{I\left( H_{i} \right)}\beta\omega_{i}}} \right)}{{Hp}_{i} = {H_{i} \cdot \left( {1 + {\Delta_{i} \cdot \left( {j + {\beta \cdot \omega_{i}}} \right)}} \right)}}$ where A′ is a vector, coordinates of which are respectively formed by: ω_(i)·[R(H _(i))·I(Ip _(i))−R(Ip _(i))·I(H _(i))], B′ is a vector, the coordinates of which are respectively formed by: R(H _(i))² +I(H _(i))²−(R(H _(i))·R(Ip _(i))+I(H _(i))·(Ip _(i))), for i ranging from 1 to N, where R(H_(i)) is a real part of the prediction H_(i), I(H_(i)) is an imaginary part of the prediction H_(i), R(Ip_(i)) is a real part of the at least one measurement Ip_(i), I(Ip_(i)) is an imaginary part of the at least one measurement Ip_(i), a fourth calculator for computing at each iteration the error E_(i) according to the equation E _(i) =Ip _(i) −Hp _(i), the first calculator having a field model for computing at each iteration at least one of a position of the at least one receiver and of an orientation of the at least one receiver by computing the prediction H_(i) corresponding to the criterion C.
 15. The magnetic locator as claimed in any one of claims 11 to 14, wherein the field model to compute the position and/or the orientation of the at least one receiver by computing the prediction H_(i) of the at least one measurement as a function of the criterion C uses a Levenberg-Marquardt minimization algorithm. 